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Abstract. 

We discuss the Buchert equations, which describe the average expansion of an 
inhomogeneous dust universe. In the limit of small perturbations, they reduce to 
the Friedmann-Robertson- Walker equations. However, when the universe is very 
inhomogeneous, the behaviour can be qualitatively different from the FRW case. In 
particular, the average expansion rate can accelerate even though the local expansion 
rate decelerates everywhere. We clarify the physical meaning of this paradoxical feature 
with a simple toy model, and demonstrate how acceleration is intimately connected 
with gravitational collapse. This provides a link to structure formation, which in turn 
has a preferred time around the era when acceleration has been observed to start. 
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The coincidence problem. Perhaps the most surprising discovery in modern cosmology 
is that the expansion of the universe has apparently accelerated in the recent past. 
The observations have usually been interpreted in the context of the homogeneous 
and isotropic Friedmann-Robertson- Walker (FRW) model. The Einstein equation then 
reduces to the FRW equations: 

3- = -4nG N (p + 3p) (1) 

(X 

3^ = 8nG N p - 3-| (2) 
p + 3-(p + p)=0, (3) 

CL 

where Gn is Newton's constant, dot denotes derivative with respect to the proper time 
t of comoving observers, a(t) is the scale factor of the universe, p(t) and p(t) are the 
energy density and pressure, respectively, and k — 0, +1, —1 gives the spatial curvature 
of the hypersurfaces of constant t. 

According to (1), the expansion of a universe which contains only ordinary matter 
with non-negative energy density and pressure decelerates, as expected on the grounds 
that gravity is attractive. The observations can then be explained only with repulsive 
gravity, either by introducing matter with negative pressure or by modifying the theory 
of gravity. Such models suffer from the coincidence problem: they do not explain why 
the acceleration has started in the recent past. If there is a dynamical explanation, it is 
presumably related to the dynamics we see in the universe. The most significant change 
at late times is the formation of large scale structure, so it seems a natural possibility 
that the observed deviation from the prediction of deceleration in the homogeneous and 
isotropic matter-dominated model would be related to the growth of inhomogeneities. 

The fitting problem. The rationale for the FRW model is that the universe appears 
to be homogeneous and isotropic when averaged over large scales. According to this 
reasoning, one first takes the average of the metric and the energy-momentum tensor, 
and then plugs these smooth quantities into the Einstein equation. However, physically 
one should first plug the inhomogeneous quantities into the Einstein equation and then 
take the average. Because the Einstein equation is non-linear, these two procedures are 
not equivalent. This leads to the fitting problem discussed by George Ellis in 1983 [1]: 
how does one find the average model which best fits the real inhomogeneous universe? 
The difference between the behaviour of the average and smooth quantities is also known 
as backreaction [2]. 

The Buchert equations. Let us look at backreaction in a universe which contains only 
zero pressure dust. Assuming that the dust is irrotational, but allowing for otherwise 
arbitrary inhomogeneity, the Einstein equation yields the following exact local equations 
[3]: 

+ \e 2 = - AitGnp - 2a 2 (4) 
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V = 8nG N p -l^ + a 2 (5) 
p + 6p = 0, (6) 

where 9(t, x) is the expansion rate of the local volume element, p(t, x) > is the energy 
density, a 2 (t,x) > is the shear and TZ(t,x) is the spatial curvature. At first sight it 
might seem that acceleration is ruled out just as in the FRW case: according to (4), 
the expansion decelerates at every point in space. However, let us look more closely. 
Averaging (4)-(6) over the hypersurface of constant t, with volume measure y^~g and 
volume V(t) = J d 3 x\/^g, we obtain the exact Buchert equations [3]: 



3- = -4nG N (p) + Q (7) 
a 

^ 2 =^G N (p)- l -(n)-\Q (8) 

d t (p)+^(p) = (9) 

Q = 2 -({6 2 )-{6) 2 )-2{a 2 ) , (10) 



where the scale factor a is defined as a (t) oc V(ty/ 3 and (A) = j d 3 xy^A/ j d 3 x^ 
is the spatial average of the quantity A. 

The Buchert equations (7)-(10) are the generalisation of the FRW equations (1)- 
(3) for an inhomogeneous dust universe. The effect of the inhomogeneity is contained 
in the new term Q. If the variance of the expansion rate and the shear are small, the 
Buchert equations reduce to the FRW equations. In fact, this is the way to derive 
the FRW equations (for dust) and quantify their domain of validity. The applicability 
of linear perturbation theory also follows straightforwardly, since Q is quadratic in 
small perturbations. Conversely, the Buchert equations show that when there are 
significant inhomogeneities, the FRW equations are not valid. In particular, the average 
acceleration (7) can be positive if the variance of the expansion rate is large. 

It seems paradoxical that the average expansion can accelerate even though the 
local expansion decelerates everywhere. We will clarify the physical meaning of the 
average acceleration with a toy model, and find that the answer sounds as paradoxical 
as the question: the average expansion rate accelerates because there are regions which 
are collapsing. 

Acceleration from collapse. When the universe becomes matter-dominated, it is well 
described by the FRW equations plus linear, growing perturbations. Within an 
overdense patch, the expansion slows down relative to the mean and eventually the 
perturbation turns around, collapses and stabilises at a finite size. Underdense patches 
in turn expand faster than the mean and become nearly empty regions called voids. Part 
of the universe is always undergoing gravitational collapse, and part is always becoming 
emptier. 

When a density perturbation becomes of order one, it passes outside the range of 
validity of the FRW equations plus linear perturbations. A simple non-linear model for 
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structure formation is the collapse of a spherical overdensity. In this case the Buchert 
equations (7)-(10) reduce in the Newtonian limit to the spherical collapse model [4], 
according to which the overdense region behaves like an independent FRW universe 
with positive spatial curvature (see e.g. [5]). Likewise, an underdense region can be 
described as a FRW universe with negative spatial curvature. 

What is true for the perturbations also applies to the mean: once there are large 
perturbations in a significant fraction of space, the average expansion no longer follows 
the FRW equations, opening the door for accelerating expansion. 

Let us consider a toy model which consists of the union of an underdense region 
representing a void and an overdense region representing a collapsing structure. We 
denote the scale factors of the regions a±, a 2 , where 1 refers to the underdense and 2 to 
the overdense region; the corresponding Hubble parameters are H\ = di/ai, H 2 = a 2 /a 2 . 
The overall scale factor is given by a 3 = a\ + a 2 , and the average Hubble and deceleration 
parameters are 

3 3 

H=-= + ^H#2 = #1 (1 - v + vh) (11) 

a a\ + a A 2 a\ + a\ 

J_a_ 1-v vh 2 v(l-v)(l-h) 2 

q ~ H 2 a~ qi (l-v + hv) 2+>l2 (l-v + hv) 2 (1-v + hv) 2 ' { } 

where qi, q 2 are the deceleration parameters of regions 1 and 2, v = a\j (af + a 3 ,) is the 
fraction of volume in region 2, and h = H 2 / Hi is the relative expansion rate. 

The Hubble rate is simply the volume- weighed average of Hi and H 2 . This is not 
the case for the deceleration parameter q: there is a third term which is negative and 
contributes to acceleration, corresponding to the fact that Q in (7) is positive. 

Let us for simplicity take the void region to be completely empty, so that a± oc t, 
qi = 0. For the overdense region we have a 2 oc 1 — cosg, where 9 is the development 
angle determined by t oc 9 — sing. This gives q 2 = (1 — cos 9)/ sin 2 6. The overdense 
region expands from 8 = 0, turns around at 9 = n and is taken (by hand) to stabilise 
at 9 = 37r/2. Denoting the volume fractions occupied by the regions at turnaround 
fi = 1 - /2, /2, we have 

/ 2 tt 3 (1 -cos#) 3 u _ sin 9(9 - sing) 

" _ M /i\ 9 • (13) 



8(1 -f 2 )(9- sin^) 3 + / 2 7r 3 (l -cosfl) 3 ' (1-cos! 

Inserting (13) into (12), it is straightforward to establish that the acceleration can 
be positive. Figure 1 shows the deceleration parameter q for this toy model with f 2 = 0.3 
and for the ACDM model with = 0.7 at 9 = 3tt/2. The toy model is not to be taken 
seriously beyond qualitative features, but we see that it is possible for backreaction to 
produce acceleration not dissimilar to what is observed. 

The explanation of the paradox that collapse induces acceleration is simple. The 
overdense region slows down the average expansion rate (11). As the volume fraction 
of the overdense region decreases, the contribution of the underdense region will 
eventually dominate, and the expansion rate will rise. This clarifies the equations 
(7), (10): larger variance of the expansion rate means that the volume fraction of the 
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Figure 1. The deceleration parameter q in the 
toy model (blue, solid) and ACDM (red, dashed). 



fastest expanding region rises more rapidly, contributing to acceleration. In fact, the 
deceleration parameter q can become arbitrarily negative: q < — 1 does not require 
violation of the null energy condition, unlike in FRW models. 

While the overdense region is essential for acceleration, it is the underdense region 
which raises Ht above the FRW value: unlike in ACDM, these two effects are distinct. 
The fact that gravity is attractive implies that acceleration driven by inhomogeneities 
satisfies Ht < 1 [6]. One consequence is that acceleration cannot be eternal, as 
is physically clear: acceleration erases inhomogeneities, so their effect decreases, 
terminating acceleration. However, inhomogeneities can then become important again. 
In the real universe, where perturbations are nested inside each other with modes 
constantly entering the horizon, this could lead to oscillations between deceleration 
and acceleration. 

Oscillating expansion could alleviate the coincidence problem, but structure 
formation does also have a preferred time near the era ~ 10 billion years where 
acceleration has been observed. For cold dark matter, the size of structures which 
are just starting to collapse, relative to the horizon size, rises as structure formation 
proceeds, saturating once all perturbations which entered the horizon during the 
radiation-dominated era have collapsed. The size of structures becomes of the order 
of the maximum size around 10-100 billion years. This is encouraging with respect to 
the coincidence problem, as one would expect the corrections to the FRW equations to 
be strongest when collapsing regions are largest. 

Outlook. The possibility that sub-horizon perturbations could lead to acceleration has 
been explored before [7], and backreaction has been demonstrated with an exact toy 
model [8]. However, the conceptual issue of how average expansion can accelerate even 
though gravity is attractive -how gravity can look like anti- gravity- had been murky. 

We have discussed how acceleration is intimately associated with gravitational 
collapse. This provides a link to structure formation, which in turn has a preferred 
time close to the observed acceleration era. Backreaction could thus provide an elegant 
solution to the coincidence problem, without new fundamental physics or parameters. 
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To see whether this promise is realised, we should quantify the impact of structure 
formation on the expansion of the universe to make sure that we are fitting the right 
equations to increasingly precise cosmological observations. 

References 

[1] Ellis G F R, Relativistic cosmology: its nature, aims and problems, 1983 The invited papers of 
the 10th international conference on general relativity and gravitation, p 215 
Ellis G F R and Stoeger W, The 'fitting problem' in cosmology, 1987 Class. Quant. Grav. 4 1697 

[2] Ellis G F R and Buchert T, The universe seen at different scales, 2005 Phys. Lett. A347 38 
[gr-qc/0506106] 

[3] Buchert T, On average properties of inhomogeneous fluids in general relativity I: dust cosmologies, 

2000 Gen. Rel. Grav. 32 105 [gr-qc/9906015] 
[4] Kerscher M, Buchert T and Futamase T, On the abundance of collapsed objects, 2001 Astrophys. 

J. 558 L79 [astro-ph/0007284] 
[5] Padmanabhan T, Structure formation in the universe, 1993 Cambridge University Press, 

Cambridge p 273 

[6] Rasanen S, Backreaction and spatial curvature in dust universes, 2006 Class. Quant. Grav. 23 

1823 [astro-ph/0504005] 
[7] Rasanen S, Dark energy from backreaction, 2004 JCAP0402(2004)003 [astro-ph/0311257] 

Rasanen S, Backreaction of linear perturbations and dark energy [astro-ph/0407317] 
[8] Rasanen S, Backreaction in the Lemaitre-Tolman-Bondi model, 2004 JCAP041 1(2004)010 [gr- 

qc/0408097] 



